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ABSTRACT 

We present a new set of dissipationless N-body simulations aiming to better un¬ 
derstand the pure dynamical aspects of the “Fundamental Plane” (FP) of elliptical 
galaxies. We have extended our previous hierarchical merger scheme by considering the 
Hernquist profile for the initial galaxy model. Two-component Hernquist galaxy mod¬ 
els were also used to study the effect of massive dark halos on the end-products char¬ 
acteristics. We have also performed new collapse simulations including initial spin. We 
found that the one-component Hernquist mergers give results similar to those found 
for the one-component King models, namely both were able to build-up small scatter 
FP-like correlations with slopes consistent with what is found for the near-infrared FP 
of nearby galaxies. The two-component models also reproduce a FP-like correlation, 
but with a significantly steeper slope. This is in agreement with what has been found 
for elliptical galaxies at higher redshift (0.1 < z < 0.6). We discuss some structural 
properties of the simulated galaxies and their ability to build-up FP-like correlations. 
We confirm that collapses generally do not follow a FP-like correlation regardless of 
the initial spin. We suggest that the evolution of gradients in the gravitational field 
of the merging galaxies may be the main ingredient dictating the final non-homology 
property of the end products. 

Key words: galaxies: elliptical - galaxies: fundamental parameters - methods: nu¬ 
merical 


1 INTRODUCTION 


luminosity ratio (M/L) of the systems (c.f. |Ca£el^Oj_^ 


The origin ot the ’Tundamental Plane" relation (hereon I P) 
of elliptical galaxies is still unknown, despite of all the efforts 


Carvalho fc Carlberg 


to understand it since it w as discovered (Djorgovski & Davii 


1987, Dressier et al. 1987). 


On fundamental grounds, the simplest version of the 
virial theorem applied to galaxies predicts that they should 
form a family of objects following a simple projected re¬ 
lation, involving structural and kinematical variables, for 
instance: re = Cvircr^le^ ■ In this equation, Cvir is a 
structural-kinematical parameter, cro is the central projected 
velocity dispersion, le, the average surface brightness within 
the effective radius in linear units, and Ve is the effective ra¬ 
dius. The coefficient Cvir relates physical (3-D) to projected 
variabl es, like the velocity dispersion and mass distributions. 


1995, hereon CdCC95; pantas et al 


2001, hereon DCdCROl). 


; ^ 


We define a family of homologous galaxies as virialized 
systems where the kinematical and structural coefficients 
are simply constant for all galaxies, or may change but in 
a constant ratio throughout the family. If furthermore M/L 
is constant for all galaxies (or equivalently, M and L may 
change among them but in a constant ratio), then Cvir is a 
constant for a given homologous family. 

The expression of the FP is similar to that expected 

from the projected virial relation but with significantly dif- 

A tB 


ferent exponents and small scatter throughout: r, 
(where the exponents are: A ^ 1.53, B ~ —0.79, e.g. 




Pahre 


Hence, Cvir = Cvir{Cr,Cv, M/L) depends on kinematical Djorgovski & de Carvalho 1998|). Thus, in the case of ellipti- 


{Cv) and structural (Cr) coefficients, as well as on the mass- cals, it is inferred that Cvir must vary monotonically among 
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these galaxies if one is willing to retrieve the virial relation _ sideration a central density peak. Recent studies (e.g. Ger- 


for these systems. 


There are several reasonable alternatives to explain why 
and how Cvir should vary in order to explain the discrep¬ 
ancy between the virial theorem and the FP. One of them 
assumes a systematic variation of M/L with the total mass 
of the system, preserving homology. This would be r e spon- 


sible for the whole C-uir variation (e.g. Djorgoys^ 


Diogorvski fc Santiagc 


1993 


Renzini et al 


1988 


1993, Pahre & 


Djorgo4~ 1997). This dependence would be caused by sys¬ 


tematic stellar population (e.g. mean stellar age, m e tallic - 


ity, etc.) variations with mass. Pahre & Djorgovski (1997) 


have shown that there is a dependence of the FP tilt with 
wavelength, namely M/L oc where P varies with the 
photometric band (A) in which the luminosity is measured. 
This means that P decreases from the B to the K band, al¬ 
though never reaching the homologous, virial expectations. 
As discussed by Pahre & Djorgovski, the trend of P with A 
cannot be explained solely by either stellar population mod¬ 
els or non-homology (see their Figure 2). They conclude 
that a more complete scenario to explain the FP tilt has 
to invoke contributions from both effects. Broken homology 
can be achieved both in dissipationless hierarchical merg¬ 
ing scenarios (e.g., CdCC95) and in dissipative mergers of 
star-forming and gas-rich spirals, wher e the roles o f star for¬ 


mation histories are emphasized (c.f. Bekki 1998). A third 


line of reasoning for explaining the FP assumes that a more 
refined formulation to describe the equilibrium condition of 
the luminous component of the elliptical galaxies is adopting 
a “two-component virial theorem”, which assumes of course 
t hat ellipticals are dynamically dominated by a dark halo 


(Dantas et al 2000). 


In the present work, we study the origin of the FP 
tilt under the assumption that elliptical galaxies are more 
closely described as non-homologous virialized systems. 


hardt et al. 


200C 


Siopis et al. 2000) have demonstrated, 


however, that the presence of a central peak (or even the 
presence of a central black hole) should be much more com¬ 
mon in elliptical galaxies than previously thought. Here we 
consider the hierarchical merge of Hernquist models, which 
present a central density “cusp”. Second, CdCC95 only con¬ 
sider one component models. However, it is important to 
understand the effect of the halo in the dynamics of merg¬ 
ing and the consequences of its influence in the equ ilibrium 
conditions of the whole system ( Dantas et al. 2000 ). In this 
paper we consider the merge of two-component Hernquist 
models up to two generations. 

One point not addressed by DCdCROl was the initial 
difference in spin parameters between proto-galaxies and 
how that can introduce non-homology into the structural 
properties of the final objects. In order to address this point, 
we investigate how the spin parameter influences the FP of 
the collapsed objects. However, the issue which we leave for 
future work is the study of two-component collapses. This is 
an important problem, since in the currently accepted cos¬ 
mological scenarios, the luminous component collapses in 
the dark matter halo already virialized some time after the 
epoch of equality of matter and radiation energy densities. 
Our present goal is to establish the behaviour of only one- 
component collapses before analysing two-component ones, 
which can be studied under a more general approach as, for 
instance, drawing the models from high resolution cosmo¬ 
logical simulations. 

This paper is organized in the following way: in Section 
2, we present the simulation setups and initial condition 
grids; in Section 3, the end products of the simulations are 
considered in the context of the FP space and the resulting 
relaxation histories. Finally, in section 4, we discuss our main 
results. 


Hjorth & Madsen] [1995|, Capelato, de Carvalho & Carl- 

berg 

|l£|971JCiotti, Lanzoni & Renzinil [19961 [Busarello et al. 

1997 

, Graham & Colless 

1997 

, Bekki [1998|). In a hierarchical 


galaxy formation scenario, galaxies are built-up by suces- 
sive merge of larger and larger systems. Recent observations 
have reinforced the idea of hierarchical merger a s a reason- 
able mecha nism to form elliptical galaxies (e.g. Bender & 
Saglia : 999|) , although dissipation seems to be an important 
ingredient. In any case, numerical investigations of merg¬ 
ing seem to be fundamental in understanding the scaling 
relations of these objects. The numerical work of CdCC95 
has shown, for instance, that the FP correlations can arise 
naturally from objects that are formed by dissipationless 
hierarchical mergers of galaxies. The end products of their 
simulations result in a non-homologous family of objects, 
being the peculiar non-homology mainly determined by the 
parameter Cv varying systematically with the initial orbital 
energy of the galactic pairs. In a subsequent investigation, 
DCdCROl have shown that one-component, equal mass col¬ 
lapses of several different initial models and collapse factors 
produce nearly homologous families of objects. This result 
led DCdCROl to suggest that the driving mechanism pro¬ 
ducing non-homology would be that of merging per se. 

We extend the previous dissipationless numerical inves¬ 
tigations in several aspects. First, the equilibrium models 
considered by CdCC95 (King spheres) do not take into con- 


2 SIMULATIONS SETUP AND DEFINITION 
OF CHARACTERISTIC PARAMETERS 

2.1 Computer Plataforms and Codes 


The simulat ions were run using tw o C translations (c.f. Du- 
binski 198^ ) of the TREECODE (Barnes & Hut 1986): a 
non-parallel version, which was used to run less CPU time 
consuming, one-component model simulations; and a paral¬ 
lel version, run for two-component model simulations. The 
computational plataforms used were: (i) For the non-parallel 
code: Workstations Sun-Sparc; Sun-Ultras (1, 2, 5, 10, and 
30); and Sun E250; (ii) For the parallel code: Silicon Graph¬ 
ics Origin 2000 with four processors using MPI (“Message 
Passing Interface”), IRIX operational system; and a “clus¬ 
ter” formed by four Pentium HI, 650 MHz, working in par¬ 
allel using LAM (“Local Area Multicomputer”) 6.3.2/MPI 
2 C-I--I-, Linux. 


Quadrupole correction terms, according to Dubinski 


were used in the force calculations for all simulations. 
In Table we list the main parameters of the simulations 
setup adopted in this work. These parameters were carefully 
chosen in order to conform to the constraints of resolution 
and collisionlessness given the total number of particles used 
in each type of simulation (more details for the choice of pa- 
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Table 1. Initial Parameters of the Codes 


Parameter 

Value 

6 : tolerance 

0.8 

e: softening 


^ Non-parallel code 

0.05 

^ Parallel code 


Luminous component 

0.07 

Dark component 

0.7 

At: time integration step 

0.025 


rameters can be found in CdCC95 and DCdCROl). In par¬ 
ticular, the choice of the number of particles was also based 
on the operational constraints due to CPU times. Merging 
generaly involves CPU time-consuming runs for it includes 
the evolution since the initial orbital phase, before the effec¬ 
tive merge of the systems. This forced us to use a relatively 
small number of particles to cover a wider grid of initial con¬ 
ditions. These numbers, however, are well above the lower 
bound discussed in DCdCROl. 


2.2 Initial Condition Grid of the Models 

2.2.1 Computational Units 

The units used in our simulations were all set to match those 
of CdCC95 and DCdCROl: the mass and lenght units were 
set to Mu = 10^° Mq and Lu = 1 kpc, respectively. These 
values, and G = 1, fix our time and velocity units to Tu = 
4.72 Myr and vu = 207 km s“^, respectively. 


2.2.2 The Merger Models 

The initial equilibrium models for the merger pair were each 
obtained from N particles random realizations of spheres in 
hydrostatic equ ilibrium, obeying the Hernquist profile (c.f. 
Hernquist 199C). We considered both one as well as two 
component models, in equilibrium in the common potential. 

The reasons for the choice of a Hernquist profile for the 
luminous and dark components were based on the desire 
to test whether models including a central density ’cusp’ (in 
this cas e, the Hernquist models provided us with this charac¬ 


teristic’ could also reproduce the results by CdCC95. Since 


the FP parameters refer to central (effective) quantities, the 
idea was to test whether the results changed sensibly or not 
with the inclusion of an initial ’cusp’ in the models. In par¬ 
ticular, the reasons for the adoption of the Hernquist profile 
also for the halo (instead of, e .g., a truncated isothermal 
profile, e.g. Walker et al. (1996)) comes from the fact that 


the density profile behaves as ~ r~ at small radii resem¬ 
bling the Navarro, Frenk & White (1997) ‘universal’ profile, 
which results from cosmological simulations. 

The Hernquist models were truncated at a specified en¬ 
ergy cut-off: 10% least bound particles of the model were 
eliminated. Hence, the original Hernquist one-component 
model had a mass of IMu. The truncated model resulted in a 
mass of 0.9Mu. The same was applied to the two-component 
models, where both the luminous as well as the dark compo¬ 
nent were truncated by the same factor (10%): the luminous 


component has a mass of 0.9Mu and the dark component, 
a mass of 9Mu. 

We assigned to the one-component Hernquist mergers 
the labels: D, E e F, according to which generation they 
belong (D: first, E: second, and F: third). The total initial 
values of mass and number of particles of the D mergers 
were: Mtot.o = 1.8 Mu and Ntot.o = 8194, respectively 
(these values refer to the sum of the two initial merging 
models, not to one model alone). 

The two-component mergers were assigned with la¬ 
bels Z (Z01-Z09: first generation; Z10-Z13: second genera¬ 
tion). We chose the initial luminous [Ml) to dark (Mh) 
mass ratio of the initial two-component Hernquist m o del a s 
p.init = Ml/Mh = 0.1 (the results of Mihos et al. (1996) 
favor Mh ~ (4 — 8) Mdiak+buige for NGC 7252, suggest¬ 
ing our mass ratio is reasonable). The total initial mass 
of the Z mergers was Mtot,z = 19.8 Mu, with a total of 
Ntot.z = 9000 particles. Each initial two-component Hern¬ 
quist model therefore has a luminous mass of 0.9 (2250 par¬ 
ticles) and a dark mass of 9 (also 2250 particles). Note that 
since the number of particles per component is the same, 
the mass per dark matter particle is greater than that of 
the luminous particle by a factor 10. 

The initial ratio of the effective (half-mass) radius of 
dark matter to that of the luminous component was an = 
10 au. Here we briefly discuss the reasoning for choosing 


these ratios. 

ies rp ~ (4 — 7)Rd, where ro is the halo core radius of the 


Salucci & Burkert (200C) find for disk galax- 


Burkert (1995) model (ro is of same order as r^, the core 
radius of the modified isothermal model). Rd is the disk 
scale radius. Noticing that the effective radius for spirals, 
(Re), is approximately related to Rd by {Re) ~ 1.2Rd, then 
ro ~ (3.3 — 5.8)(7?f). Noticing also that (i?f) ~ Rf for 
L — Lt galaxies, where Rf is the effective r adius for giant e l- 
lipticals, and that the Re ~ 1.8153 au (c.f. Hernquist 1990 ), 
where ol is the scale radius of the Hernquist profile for the 
luminous component, one can infer that the results of Salucci 
& Burkert imply ro « (6— 10)a_L. Assuming ao is of the same 
order as oh, the scale radius of the Hernquist profile for the 
dark matter component, there is a compatibility between 
our adopted values for the initial ratio {an = 10 au) and 
the results by Salucci & Burkert (although thei r analysis was 
based on spiral galaxies). Serhardt et al, (2001), on the other 
hand, find that re,h ~ 1-2 Re for EO elliptical s, where re.h 
is the “min imum halo model” core radius (c.f. Kronawitter 
et al. 2000). Again, this can be translated to re,h ~ 2.2 au. 
It is not at all clear the correspondence between rc,/i and 
Oh, but if they have the same order of magnitude, it would 
seem to imply our value {oh = 10 au) would be somewhat 
higher than adequate. On the other hand, however, there 
are some works on the morphology and kinematics of tidal 
tails of merger models, where some inferences can be made 
on the halo prope rties by a comparison with simulations. 
Mihos et al. (1996), for instance, study models with ratios 


of mass and radius within the range of our model. They find 
a good fit to NGC 7252, favoring relatively compact, low- 
mass halos for the progenitors of the merger. Although their 
results are somewhat idealistic, our models do not seem to 
be imcompatible with what is usually adopted in the liter¬ 
ature. However, in face of the uncertainties for a reasonable 
value for the effective (half-mass) radius of dark matter to 
that of the luminous component, we check the dependency 
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of the FP results on the choice of this ratio. To that end, 
have run two sets of nine simulations similar to the Z mod¬ 
els, but using a more compact halo, namely: an = 2> ah 
and an = 5 an. These models are labeled Z01b-Z09b and 
Z01c-Z09c, respectively. 

The initial merging conditions were characterized ac- 
cording to a genera l izatio n of a prescription described in 
Binney & Tremaine (1987) (c.f. CdCC95). In this formula¬ 
tion, the initial orbit of the binary galaxy system is charac¬ 
terized, essentially, by the energy and angular momentum 
of the Keplerian orbit of two point masses equivalent to the 
initial galaxies. We defined the dimensionless energy and 
angular momentum of the orbit as: 


P _ Eorb 

( 1 ) 


Y _ ^orb 

L — ^ 

( 2 ) 




with {v^) = ^{v\){v\), Vh = .yrhirh 2 , where rui (i = 1 , 2 ) 
is the half-mass radius of the system i, and fi is the reduced 
mass of the system. A third parameter depends only on the 
dynamical structure of the initial galaxies: 


2GM 

rh{v^) 


( 3 ) 


which presents a not very large variation ( $ 20 %) among 
the initial models (A ~ 17). 

The initial separation of the models was chosen as ~ irh 
for the parabolic and hyperbolic orbits, and the apocenter 
position for the closed orbits. These initial separations were 
chosen considering that they should not be too close (imply¬ 
ing that tidal effects would be artificially disregarded due to 
the spherical symmetry of the initial models) nor too far 
away, so that time consuming CPU runs were avoided. 

By using this grid of initial conditions, the models 
merged and evolved up to ~ 30 “crossing times” (Ter = 
(jjy^ 5 / 2 /( 2 |j^|) 3 / 2 ), quantities like half-mass radius 

(r/,) and virial ratio (/3^ = ■^^) indicated no significant 
variation of the resulting system {Arn/rh 0.5%, and 
A/3„//l„ < 1% after ~ 10 Ter). 

In Tables (one-component models) and ^ (two- 
component models) we list the initial condition grids of the 
merger simulations. The two-component merger simulations 
using a more compact halo than the Z models (i.e., the Zb 
and Zc models, were an = 3 ul and an = 5 ul, respec¬ 
tively) are also included in Table 1^. We also list in Table 
^ the simulations performed by CdCC95, including several 
third generation simulations not previously published (the 
total number of particles for the first generation of these 
mergers is 8192). 


2.2.3 The Collapse Models 

In Table ^ we include the simulations performed by 
DCdCROl for easy reference. Details of the collapse models 
can be found in DCdCROl. Here we give a brief summary 
of these collapse simulations: three different initial collapse 
models were used (labeled K, A and C). All the initial mod¬ 
els have total mass M = 2Q Mu and radius R = 20 Lu, 
except the C models, which have R = 100 Lu. The K 


models were constructed from 8192 Monte Carlo realiza¬ 
tions of a spherical isotropic King model. The A models 
were constructed from sphe rical models of1638 4 Monte 
Carlo particle realizations ( Aguilar fc Merritt 199C). The C 
models were constructed according to [Carpintero fc Muzzio 
( 1995| ), with 4096 particles. The initial velocities of these 
models had gaussian profiles. All models were pertubed ac¬ 
cording to the collapse factor parameter /3 (0 /3 1, 

where /3 = 2A'o/|Ho|; Ko is the initial kinetic energy and 
Ho, the initial potential energy of the system). To the C 
models, a Hubble flow assuming Hq = 65 km s“^ Mpe”*^ 
has been incorporated. We generically denoted “cold” col¬ 
lapses those resulted from /3 —> 0, and “hot” collapses those 
resulted from /3 —> 1. 

We have included also two sets of collapse simulations 
with a range of initial solid body rotation, not discussed 
in DCdCROl. We have included spins to the unperturbed, 
initial A model in order to study their effects in the final 
systems. The reason to focus on the A models is because 
these collapses spread in the FP space, contrary to the K 
models. Although the C models are more ‘realistic’ (they 
evolve from small pertubations in a Hubble expansion), be¬ 
cause of being dinamically more complex we have avoided 
them our analysis of the spin effect on the FP (see details 
in DCdCROl). 

The method we assum e here is inspired on that of 
Wilkinson fc James| (1982). We have given a solid body 
angular velocity, ui, to each particle of the unperturbed 
A model. The value of to was chosen such that the re¬ 
sulting total kinetic energy after including the spin was 
a fraction 7 greater than the initial total kinetic energy 
(i.e., without the rotational motion). In other words, 7 = 
\Kprog — Kprog,spin\/Kprog. For the first Set of collapse simu¬ 
lations with spin, which we label ASl models, the rotational 
perturbation chosen was small, 7 ~ 5%. The total veloc¬ 
ity squared of each particle was then reduced by a range 
of P factors, producing 9 spin models with different col¬ 
lapse fators. These collapses can be directly compared to 
the A collapses studied by DCdCROl. For the second set of 
models, 7 was chosen in order to impose a maximum per¬ 
turbation to the A progenitor such that, after reducing the 
velocity field by the “hottest” perturbation we are consider¬ 
ing (viz. log P = — 0 . 01 ), the resulting model was barely able 
to collapse (total binding energy was ~ —0.007). The value 
the perturbation in this case was 7 ~ 38%. The perturbed 
progenitor was “cooled” by the same P factors as the ASl 
models. This second set of models was labeled AS3 models. 
These new collapse simulations are listed in Table 

Note that the structure (viz. potential energy) of the 
A models used here to construct the spin models did not 
allow the inclusion of a higher initial spin than that of the 
AS3-09 (log/3 = —0.01) model without disrupting the sys¬ 
tem (viz., expanding it instead of making it collapse). Higher 
spin rates could have been used, but that would imply chang¬ 
ing the structure of the progenitor by, e.g., reconfiguring 
the positions of the particles (viz., by decreasing the grav¬ 
itational radius of the system) or incrasing the total mass. 
That would change considerably the structure of the pro¬ 
genitor and would not allow a direct comparison with the 
collapse models of DCdCROl. Hence, these collapses with 
spin just represent models were an initial rotational “pertur¬ 
bation” was applied to the system. This alowed us to keep 
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Figure 1. Results of the merger simulations in terms of the FP 
paremeters. In the case of two-component models, only the lumi¬ 
nous component is presented. The stellar symbol represents the 
progenitor. The continuous line represent the prediction of the 
virial theorem for homologous systems. 


the same initial structure of the progenitor of the A mod¬ 
els, used by DCdCROl, and still make the resulting model 
collapse according to the (3 factor. 


3 THE FUNDAMENTAL PLANE OF END 

PRODUCTS 

3.1 The FP space 

We follow the method given by CdCC95 to compute the 
characteristic FP variables (re, ^e = —2.51og7e and uq) 
of the simulated models. The variables no and fie were 
combined in the vertical axis according to the usual rep¬ 
resentation of the FP projected onto the cartesian plane 
log re X log (To + ft fie- lu all the cases the 3-variate best-fit 
solutions for a plane gave (3 = 0.2 to within 10%, so we de¬ 
cided to keep this coefficient fixed at 0.2 in order to find the 
orthogonal least squares solutions for the other coefficients, 
viz. the slope a (the FP ‘tilt’), and the intercept of the fitted 
plane. 

Before analysing the final simulated models in the FP 
space, however, it is worth to comment that they reproduce 
the general structural characteristics of elliptical galaxies, 
e.g. projected triaxialities (from EO to E5 elliptical objects) 
and surface density profiles (following the Sersic law). A 
detailed discussion on t he structural p roperties of the simu- 

(|200l|). 
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lated models is given in 

In Fig. ^ we present the characteristic FP parameters 
of the objects resulting from the merging of one (D, E e 
F mergers) and two (Z mergers) component models. In the 
case of two-component models, the data shown in this fig¬ 
ure are relative to the luminous component. The best fit 
values of the FP slope (q) found for these simulations are 
indicated in the figures. The continuous line (a = 2) rep¬ 


resent the prediction of the virial theorem for homologous, 
constant M/L, systems. In Table we present the results of 
the best fit values of the data here discussed as well as the 
results obtained by CdCC95 and DCRdCOl. The results in¬ 
dicate that one-component Hernquist mergers (D, E, F) also 
reproduce reasonably the FP tilt of the elliptical galaxies, 
consistent with the results obtained with the King models of 
CdCC95. That is, for both cases, the FP slopes are consis¬ 
tent, within the errors, to that observed for infrared FP of 
nearby galaxies. tha. t is, a = 1.53 ± 0.08 (Pahre, Djorgovski 


fc de Carvalho 1998 ). 

On the other hand, the luminous/barionic objects re¬ 
sulting from the the two-component mergers form a fam¬ 
ily with a steeper relation in comparison to one-component 
mergers. (It is interesting to note that if we consider an 
equivalent to the FP space but for the dark halos of these 
merger remnants, we find that they constitute an approx¬ 
imate homologous family of objects, as indicates the value 
a = 1.872, in Table ^) In order to test the effects of a 
more compact halo on these results, we ran two groups of 
two-component Hernquist mergers with different ratios for 
the halo to luminous radius (as discussed in Section 2.2.2). 
Unlike the Z models, these runs were not followed for sub¬ 
sequent (viz., second, etc) generations because of CPU time 
limitations. In Fig. we show how these more compact halo 
mergers distribute in the FP space. The arrows over the dot¬ 
ted lines in that figure represent the range occupied by the 
first generation Z models (arr/aL = 10), for comparison. 
It is interesting to notice that the luminous component of 
the most compact halo models (Zb models), with most neg¬ 
ative Eorb’s, tend to cluster in the FP space in a similar 
manner as the K collapses. All other models tend to spread 
out sensibly along a FP-like relation. Also, the luminous 
body of these models (Zb and ZC models) tend to settle 
into systematically lower values of Ve and at higher values 
of y (= log (To + 0.2{fie)) than the Z models (an/cLL ~ 10). 
The FP ‘tilts’ of these models suggest a marginally steeper 
‘tilt’ than the Z models. 

As already mentioned, we have performed two groups 
(ASl and AS3 models) of collapse simulations with initial 
spin in order to verify the effects of the inclusion of rota¬ 
tion on the results by DCdCROl, where evidences for ho¬ 
mology were found for pure collapses. The resulting FP 
‘tilts’ for both groups (qasi = 2.204 ± 0.158; oass = 
2.306±0.250) suggests that the resulting models are slightly 
non-homologous, but in the opposite (viz. a > 2) sense from 
the observed FP ‘tilt’ of elliptical galaxies (c.f. Fig. ^). 

All these new collapse models evolved for more than 2 
Gyr (~ OOTcr), however, the “hottest” models (viz., ASl- 
09 and AS3-09, both with initial log2A/|IU| = —0.01) still 
presented a virial ratio oscillating around 2KIW ~ —1.4 
by that time. These “hottest” models seem to evolve very 
slowly and still did not reach complete virial equilibrium 
after 2 Gyr, whereas all the other models were already well 
virialized. Removing these “hottest” collapse models results 
in QiAsi— 0.01 = 1.966±0.270 and oass-o.oi = 2.190±0.349, 
which are compatible with homology. 

Although we cover a reasonable range of /3’s, it is not 
clear whether the inclusion of more intermediate /3 collapses 
would necessarily improve the statistics (i.e., decrease the 
error bar of the fit), since the “coldest” models tend to clus¬ 
ter in the FP space. On the other hand, the inclusion of 
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Table 2. Mergers of Hernquist Models (one-component) 


!st. Generation 2nd. Generation 3rd. Generation 


Run 

E 

L 

Sep. 

^part 

Run 

E 

L 

Sep. 

^ part 

Progen. 

Run 

E 

L 

Sep. 

^ part 

Progen. 

DOl 

0.0 

0 

15.0 

8194 

El 

0.0 

0 

16.0 

15777 

D01-D02 

FOl 

0.0 

0 

22.0 

31260 

E02-E02 

D02 

-3.0 

1 

11.2 

8194 

E2 

-2.0 

1 

46.4 

16064 

D03-D04 

F02 

-1.0 

1 

25.0 

30209 

E02-E05 

DOS 

-1.0 

1 

33.8 

8194 

E3 

-1.0 

1 

89.1 

16248 

D07-D08 

F03 

-10.0 

1 

15.2 

28497 

E02-E04 

D04 

-7.5 

2 

4.0 

8194 

E4 

0.0 

2 

12.0 

15416 

DOl-DOl 

F04 

-7.0 

2 

21.8 

27446 

E04-E05 

D05 

-1.0 

2 

33.4 

8194 

E5 

-3.0 

0 

12.0 

15766 

D03-D03 







D06 

0.5 

2 

15.0 

8194 

E6 

-0.3 

0 

100.0 

15523 

D05-D06 







D07 

-6.9 

3 

3.4 

8194 













DOS 

-2.8 

3 

10.9 

8194 













D09 

0.0 

3 

15.0 

8194 













DIO 

-5.0 

0 

10.0 

8194 














Table 3. Mergers of Hernquist Models (two-components) 



1st 

Generation 





2nd. General 

ion 


Run 

E 

L 

Sep. 

^part 

Run 

E 

L 

Sep. 

^part 

Progen. 

ZOl 

-4.0 

0 

70.0 

9000 

ZIO 

-2 

1 

207.66 

17807 

Z06-Z07 

Z02 

-4.0 

1 

70.0 

9000 

Zll 

-1 

1 

416.04 

17507 

Z09-Z02 

Z03 

-3.0 

0 

70.0 

9000 

Z12 

-1 

1 

416.04 

17944 

ZOl-ZOl 

Z04 

-3.0 

1 

138.2 

9000 

Z13 

0 

0 

70.00 

17954 

Z01-Z04 

Z05 

-2.0 

0 

70.0 

9000 







Z06 

-2.0 

1 

207.6 

9000 







Z07 

-2.0 

2 

205.4 

9000 







Z08 

0.0 

0 

70.0 

9000 







Z09 

0.5 

0 

70.0 

9000 








ZOlb 

-4.0 

0 

70.0 

9000 

Z02b 

-4.0 

1 

70.0 

9000 

Z03b 

-3.0 

0 

70.0 

9000 

Z04b 

-3.0 

1 

138.2 

9000 

Z05b 

-2.0 

0 

70.0 

9000 

Z06b 

-2.0 

1 

207.6 

9000 

Z07b 

-2.0 

2 

205.4 

9000 

Z08b 

0.0 

0 

70.0 

9000 

Z09b 

0.5 

0 

70.0 

9000 

zoic 

-4.0 

0 

70.0 

9000 

Z02c 

-4.0 

1 

70.0 

9000 

Z03c 

-3.0 

0 

70.0 

9000 

Z04c 

-3.0 

1 

138.2 

9000 

Z05c 

-2.0 

0 

70.0 

9000 

Z06c 

-2.0 

1 

207.6 

9000 

Z07c 

-2.0 

2 

205.4 

9000 

Z08c 

0.0 

0 

70.0 

9000 

Z09c 

0.5 

0 

70.0 

9000 


Table 4. Mergers of King Models (one-component, CdCC95 plus 3rd. generation new data) 


1st. Generation 2nd. Generation 3rd. Generation 


Run 

E 

L 

Run 

E 

L 

Progen. 

Run 

E 

L 

Progen. 

R1 

0.0 

0 

HI 

0.5 

3 

R17-R17 

H14 

-2.0 

3 

R6-H3 

R2 

-4.0 

1 

H2 

-2.0 

1 

R6-R6 

H15 

-2.0 

4 

H1-H3 

R3 

-3.0 

1 

H3 

-4.0 

1 

R17-R17 

H16 

0.5 

4 

Hl-Hl 

R4 

-2.0 

1 

H4 

-2.0 

1 

R8-R8 

H17 

-1.77 

4 

H13-H13 

R5 

-1.0 

1 

H5 

-2.0 

1 

R14-R14 

H18 

-3.0 

3 

H13-H13 

R6 

0.5 

1 

H6 

-2.0 

3 

R2-R2 

H19 

-0.5 

3 

H10-H19 

R7 

-7.5 

2 

H7 

-2.0 

-3 

R2-R2 





R8 

-5.7 

2 

H8 

-2.0 

3 

R9-R9 





R9 

-1.0 

2 

H9 

-3.0 

2 

R9-R9 





RIO 

0.0 

2 

HIO 

-3.0 

-2 

R9-R9 





Rll 

0.5 

2 

Hll 

-2.0 

1 

RIO-RIO 





R12 

-7.9 

3 

H12 

-2.0 

1 

Rl-Rl 





R13 

-6.9 

3 

H13 

0.5 

2 

Rll-Rll 






R14 -5.1 3 

R15 -2.8 3 

R16 -1.0 3 

R17 0.0 3 





“Fundamental-Plane”-Like Relations 7 


Table 5. Collapses (one-component, DCdCROl) 


K Models 

A Models 



C Models 



^part 

^ 8192 

^part 

- 16384 



^part 

- 4096 







n 

= 1 

71 

- 0 

71 = 

2 

log/3 

Run 

log/3 

Run 

log/3 

Run 

log/3 

Run 

log/3 

Run 

-4.00 

KOI 

-4.00 

AOl 

-3.75 

COl 

-3.75 

Cll 

-3.75 

C21 

-3.75 

K02 

-3.50 

A02 

-3.50 

C02 

-3.50 

C12 

-3.50 

C22 

-3.50 

K03 

-3.00 

A03 

-3.25 

C03 

-3.25 

C13 

3.25 

C23 

-3.25 

K04 

-2.50 

A04 

-3.00 

C04 

-3.00 

C14 

-3.00 

C24 

-3.00 

K05 

- 2.00 

A05 

-2.50 

C05 

-2.50 

C15 

-2.50 

C25 

-2.75 

K06 

-1.50 

A06 

- 2.00 

C06 

- 2.00 

C16 

- 2.00 

C26 

-2.50 

K07 

-1.25 

A07 

-1.50 

C07 

-1.50 

C17 

-1.50 

C27 

-2.25 

K08 

- 1.00 

AOS 

- 1.00 

COS 

- 1.00 

C18 

- 1.00 

C28 

- 2.00 

K09 

-0.75 

A09 

-0.90 

C09 

-0.90 

C19 

-0.90 

C29 

-1.75 

KIO 

-0.50 

AlO 

-0.80 

CIO 

-0.80 

C20 

-0.80 

C30 

-1.50 

Kll 

-4.10 

AOlb 

-4.00 

COlb 





-1.25 

K12 

-3.60 

A02b 

-3.60 

C02b 





- 1.00 

K13 

-3.40 

A03b 

-3.40 

C03b 





-0.75 

K14 

-3.10 

A04b 

-3.10 

C04b 





-0.50 

K15 

-2.75 

A05b 

-2.25 

C06b 





-0.25 

K16 

-2.25 

A06b 

-1.75 

C07b 





- 0.01 

K17 

-1.75 

A07b 

-1.25 

C08b 







-1.25 

A08b 

-0.25 

C09b 







-0.95 

A09b 

- 0.10 

ClOb 







-0.85 

AlOb 









-0.25 

All 









- 0.10 

A12 








Table 6. Collapses (the “A” Model Progenitor with Initial Spin) 


ASl Models 

Npart ^ 16384, 7 - 5% 

AS3 Models 

Npart = 16384, 7 - 38% 

log/3 

Run 

log/3 

Run 

-4.00 

ASl-01 

-4.00 

AS3-01 

-3.50 

ASl-02 

-3.50 

AS3-02 

-3.00 

ASl-03 

-3.00 

AS3-03 

-2.50 

ASl-04 

-2.50 

AS3-04 

- 2.00 

ASl-05 

- 2.00 

AS3-05 

-1.50 

ASl-06 

-1.50 

AS3-06 

- 1.00 

ASl-07 

- 1.00 

AS3-07 

-0.50 

ASl-08 

-0.50 

AS3-08 

- 0.01 

ASl-09 

- 0.01 

AS3-09 


“hotter” collapses would only exacerbate the observed “in¬ 
verted” (viz. a > 2) non-homology. These results seems to 
indicate that an initial spin is not sufficient to produce non¬ 
homology, at least of the same nature of mergers. 


3.2 Spin Analysis 


In this section we briefly analyse the how the final spin of the 
models depend on the initial condition. We paramet rize the 
spin by the dimensioness quantity A, defined by (c.f. Peeble: 


1971| ): 

L\^ 

GM^/^ ’ 


(4) 


where L is the total angular momentum of the system about 
its baricenter, E the total energy of the system, and M the 
total mass (as already mentioned, G = 1). 

Fig. ^ shows how the spin of the mergers distribute as 
a function of the initial orbital angular momentum of the 
pre-merger pair. First, it can be seen that indeed there is a 
transfer of Lorb to the final spin of the merger, since higher 


Lords produce systematically higher final spins. Second, in¬ 
termediate Lords (1 < Lorb < 3) producc objects with 
spins compatible with boxy ellipticals. We note, however, 
that the position of the merger products on the FP depends 
very little on Lorb (c.f. CDCC95). In other words, Lorb = 3 
mergers could perfectly be produced from Lorb ~ 0 mergers, 
and the final products would have approximately the same 
positions on the FP. 

Fig. I plot both mergers and collapses as a function 
of the initial conditions Eorb and f3, respectively. It can be 
seen that mergers from a wide range of Eords are able to 
produce objects in the observed range of ellipticals, as op¬ 
posed to collapses, which fail in this respect. It is interesting 
to notice that “colder” collapses reach a higher degree of 
final spin than the “hottest” ones. This seem to imply that 
the initial rotational pertubations are amplified in the “cold¬ 
est” collapses. Yet, as we have seen, these “colder” objects 
still manage to become approximately homol^ous (see fits 
for ASl — 09 and ASS — 09 sequences in Table Q). Evidently, 
these results must be interpreted with caution, since we did 
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Figure 2. Results for the luminous component of the Hernquist 
merger simulations, using more compact halos. Left panel: Zb 
mergers Right panel: Zc mergers. The three lines on these panels 
are reproductions of the fits shown in Fig.l, for comparison. In 
particular, the solid line is the prediction of the virial theorem for 
homologous systems. The arrows over the dotted line represent 
the range occupied by the first generation Z models (a^/a^ = 
10), for comparison. 


CM 

d 

+ 


bO 

O 



Figure 3. Results of the collapse simulations with initial spin 
in terms of the FP paremeters. Left panel: ASl collapses, con¬ 
structed from the progenitor of the “A” models. This progenitor 
received a solid body rotation resulting on a ~ 5% perturbation to 
the initial total kinetic energy of the system, and then “cooled” 
by a range of (5 factors. Right panel: ASS models, produced in 
similar manner as the ASl collapses, but receiving a larger initial 
rotational perturbation (~ 38%). The three lines on these panels 
are reproductions of the fits shown in Fig.l, for comparison. In 
particular, the solid line is the prediction of the virial theorem for 
homologous systems. 


not reconfigure the initial structure of the progenitor in or¬ 
der to include higher initial spins. 


3.3 The Virial Coefficients 

We use another diagnostics for testing homology on the fi¬ 
nal simulated objects. A good quantitative measure in this 
case is the direct computation of the kinematical-structural 


Table 7. PF Best Fit Values 


Model 

a zb 5a 

^fit 

One-com 

ponent Models: 


D, E, F Mergers 

a = 1.414 ± 0.1S2 

20 

King (CdCC95) Mergers 

a - 1.36 zb 0.08 

17 

K Collapses 

no fit: cluster of data points 

17 

C Collapses: 

n = 0 

c. = 2.070 ± 0.123 

10 

n = 1 

a - 2.161 ± 0.087 

19 

n = 2 

a = 2.033 ± 0.342 

10 

A Collapses 

a = 1.954 ± 0.123 

22 

ASl Collapses (all) 

a - 2.306 ± 0.250 

9 

ASl Collapses (removing ASl-09) 

a = 2.204 ± 0.158 

8 

ASS Collapses (all) 

a - 2.190 ± 0.349 

9 

ASS Collapses (removing ASS-09) 

a - 1.966 ± 0.270 

8 

Two-com 

ponent Models: 


Z models (a// = 10 aiL)'- 

luminous comp. 

Q = 1.022 ± 0.046 

13 

dark comp. 

a = 1.872 ± 0.152 

13 

both components 

a - 1.176 ± 0.070 

13 

Zb models (a/f = S cliL)'- 

luminous comp. 

a = 1.004 ± 0.123 

9 

Zc models (ajf = 5 

luminous comp. 

c. = 1.017 ± 0.105 

9 
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Figure 4. Distribution of the spin of the mergers as a function 
of the initial orbital angular momentum of the pre-merger pair. 


or virial coefficients {Cr,Cv), as described in CdCC95 and 


DCdCROl: 


0 

III 

(5) 

and 



(6) 


where is the effective radius (the radius that defines a 
sphere containing half of the total luminosity of the sys¬ 
tem): L{< Te) = Ltotl’i- uo is the central projected velocity 
dispersion, and ra is the gravitational radius, defined by 
ra = GM^/|fF|, where W is the total potential energy of 
the system. 7e = L(< r^)l'Kr\ is the mean surface bright- 










Figure 5. Distribution of the spin of collapses and mergers as 
a function of the initial collapse factor (/3) and orbital energy 
{Earb)y respectively. 


ness within r^, in linear units. Then, 7e = Ci J , with 

Ci = . Inserting the equations above into the virial 

relation ((u^) = GM/rc), we find that re = CvirO’^I^^, 
where; 


C. 


2tvGCi ■ 


( 7 ) 


Since, by construction, Ci (viz. M/L) is constant among 
the models, the computation of Gr and Gv directly gives 
the measure of non-homology among the simulated models. 
Note that for two-component systems, re and {v^) are cal¬ 
culated from, respectively, the total potential and kinetic en¬ 
ergy of the system. Values of (Jq and correspond, however, 
only to the visible/barionic matter. As already pointed out, 
non-homologous objects are those which the kinematical- 
structural coefficients assume different values for each ob¬ 
ject. The results are presented in Fig. ^ where we plot the 
coefficients as a function of the initial conditions. 

First, we find that the structural coefficients, Cr, attain 
different ranges of values for one and two-component models: 
for one-component mergers, 2.5 ^ Gr ^ 3.5; whereas for two- 
component mergers, 8 ^ Cr $ 15. This difference is due to 
the presence of the massive halo in the two-component mod¬ 
els, which pushes the gravitational radius to larger values, 
as compared to the one-component systems. This increase of 
re cannot, however, be compensated by r^, which depends 
only on the structure of the luminous core. The kinematical 
coefficients, Cv, on the other hand, show similar ranges for 
both types of mergers. The product GrGv (c.f. upper panel 
of Fig. ® therefore attain larger values for two-component 
models than for one-component ones. 

A more relevant aspect of Fig. ^ is the fact that the kine- 
matical/structural coefficients vary in a systematic manner 
as a function of the initial orbital energy of the merging mod¬ 
els, which is in agreement with the results found by CdCC95. 
This behaviour seem to be an important feature distinguish¬ 
ing mergers from collapses. Indeed, collapses as a whole are 
approximately homologous objects, although some distinc¬ 
tions between “cold” and “hot” collapses are found (detailed 
discussion for collapses can be found in DCdCROl). There 
seems to be no correlation with the orbital angular momen- 
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One-Component Mergers 



Figure 6. Ratio of the virial coefficients Cvir, Cr and Cv as 
a function of the initial conditions. Symbols are: 5-point stars 
for the D models; 3-point stars for the E models; 7-point stars 
for the E models; 6-point asterisks for the Z models (1st. gen.); 
4-point asterisks the for Z models (2nd. gen.); open circles for 
the King models (1st. gen.); open squares for the King models 
(2nd. gen.); and open triangles for the King models (3rd. gen.). 
A dashed horizontal line is indicated and represents an arbitrary 
homologous family of objects for comparison. 


turn, as can be seen from an inspection of Fig. ^ On the 
other hand, it can be seen that the deviation from homol¬ 
ogy is more accentuated for two-component mergers: If we 
take the total fractional difference of y{Eorb) = {GrGv)/2 'k, 

— |1 /(A'o7 ‘& max) y^Eorb min) |/t/(T/or'6 max) ^ tve find dy ~ 

0.9 for two-component mergers whereas Sy ~ 0.5 for one- 
component models {Sy is ~ 0 for homologous objects). This 
quantity therefore reproduces the deviation from homology 
as pictured in the FP space (c.f. Fig. |^, with the advan¬ 
tage that it is possible to trace the source of non-homology 
from the corresponding fractional differences of the Cr and 
the Cv coefficients separately. For one-component mergers, 
S{Cr) ~ 0.14, S{Cv) ~ 0.25; for two-component mergers, 
5{Cr) ~ 0.38, 5{Cv) ~ 0.42. Therefore, for one-component 
mergers, Cv contributes more to the non-homology than Cr, 
a feature that can be seen clearly from an inspection of Fig. 

i 

3.4 The Ratio of “Central” to “Envelope” Kinetic 
Energies 

The results of the previous section demonstrate that 
the (central) non-homology effect which characterizes our 
merger simulations has a predominant kinematical origin. 
Now we will analyse the behavior of the total kinetic en¬ 
ergy interior to a given radius as compared to the corre- 
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King Mergers Hernquist One-Component Mergers 
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Figure 7. Ratio of the kinectic energies exterior and interior to 
three distinct radii: 0.2rf^, 0.5r^, e Irj, [notation of the figure: 
Kx = Ktot(r > IKtot{r < xrh), normalized by the corre¬ 
sponding value of the progenitor (unperturbed model)], as func¬ 
tion of the initial conditions. The dashed line indicates that the 
value of Kx = 1. The symbols are: diamonds for the first gener¬ 
ation, circles for the second generation, and crosses for the third 
generation. 


sponding kinetic energies exterior to that radius. In other 
words, if we call Ktot{< r) the “central” kinetic energy of 
the the system and Ktot{> r) the kinetic energy of its “en¬ 
velope”, then a measure the ratio of these quantities,= 
Ktot{r > xrh)IKtot[r < xr^), normalized to its progeni¬ 
tor value, should reveal, at least in a gross sense, the ef¬ 
fects of the process of relaxation. This process will therefore 
be viewed as alterations of the kinetic energies of the more 
gravitationally bounded (“central”) particles against the less 
bounded ones (“envelope”). Thus if Kx = 1 then the end 
product model presents the stratification of kinetic energies 
similar to the progenitor model. If however Kx < 1 then 
the “central” particles are “hotter” than the “envelope”, as 
compared to the progenitor. 

We analysed the kinetic energy ratio, Kx, as a function 
of the initial conditions (collapse factor or orbital energies, 
for the mergers), for three different radii {x = 0.2, 0.5 and 
1). The results are shown in Fig. ^ for the mergers models 
and in Fig. ^ for the collapse models. 

Most collapse models present Kx < 1 for any x. The 
“hottest” K collapses on the other hand approach Kx —> 1. 
In other words, the values of Kx do not change with the ini¬ 
tial collapse factor (/3), except for the “hottest” K collapses. 
Moreover the values of Kx are similar among different mod¬ 
els, except, possibly, the C models which seems more noisy 
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Figure 8. The same as previous figure, but for collapses. 


than the others. For larger x, all the collapse models have 
Kx 1. The general trend is that the collapse models are 
centrally “hotter” than the corresponding progenitor, inde¬ 
pendently of the initial j3 (except for /3’s very close to 1) and 
the initial model used. 

The stratification of kinetic energies in the case of merg¬ 
ers is not similar to the collapses. For mergers, it is clear 
that Kx is a systematic function of the initial orbital energy 
of the pairs. In other words, for mergers with more nega¬ 
tive initial orbital energies Kx —> 1, showing no difference 
with their progenitors, whereas the ones with less negative 
energies deviate more from the progenitor, and in a sys¬ 
tematic way, towards Kx < 1. The magnitude of the devia¬ 
tion from Kx = 1 also depends on the merging models: for 
instance, is greater for the King models, intermediate for 
the Hernquist one-component models, and smaller for the 
Hernquist two-component models. It also seems to slightly 
increase for increasingly merger generations. There are also 
examples where Kx > 1 (some Hernquist two-component 
models, with x = 0.2). In other words, the behavior of Kx 
for mergers seems to be more complex than collapses and 
shows a clear systematic dependency on the initial orbital 
energy of the pairs, in the same sense that the virial coeffi¬ 
cients depend systematically on Eorb- 

In case of mergers, the systematic dependency of Kx on 
Eorb begins to flatten and tend to be erased for sufficiently 
large values of x. This in fact shows that the merger models 
tend to a similar stratification of the “central” and “enve¬ 
lope” kinetic energies at sufficiently large radius. In other 
words, the different Kx values among the merger models 
are not only a function of the initial orbital energy but is a 
function of x as well, so that the correlation Kx x Eorb is 
stronger at the very center of the models and tend to dis¬ 
appear at sufficiently large radii. This shows therefore that 
the effect is intimately related to the central parts of the 
system. 

Our detailed description of the ratio of kinetic energies 
behaviour among models, as given in this section, seem to re¬ 
inforce the idea that the non-homolgy in mergers is a central 
effect ruled by how the particles gain kinetic energy during 
the merger. In other words, the non-homology seem to have 
a dynamical origin which is not present in simple collapses. 

In the following, we focus on the analysis of the relax¬ 
ation history of both mergers and colapses, which may help 
us to find clues for understanding the dynamical processes 
that are at the origin of the non-homology of mergers. 

















































Figure 9. The behaviour of the virial ratio 2KIW during the 
evolution of different and representative types of models. Upper 
panel: collapses. Lower panel: mergers. 


3.5 Relaxation Histories 

3.5.1 Evolution of the Virial Ratio 2KfW 

In order to trace a measure of the fluctuations of the gravita¬ 
tional potential on its way to equilibrium, we compared the 
behaviour of the virial ratio 2KfW (measured for the whole 
system, including escapers) during the evolution of differ¬ 
ent and representative types of models, namely: a “cold” 
(AOl model, log/3 = —4) and a “hot” (A09 model, log/3 = 
—0.75) collapse; against a “rapid” (DIO model, Eorb = —5, 
Lorb = 0) and a ’’slow” (DO model, Eorb = 0, Lot-6 = 3) one- 
component merger. Fig. H shows the behaviour of these rep¬ 
resentative models. Notice that some models do not stabilize 
around 2KIW = —1, as would be expected for a virialized 
model. This is due to the fact that we are measuring the 
virial coefficient using the complete particle data, including 
particles with positive energies which have escaped the sys¬ 
tem. More “violent” relaxations produce more escapers, and 
the resulting virial ratio stabilizes around some other value 
slightly different than —1. 

We notice that the most rapidly merging system suffer 
only one major fluctuation of 2K/W, subsequently rapidly 
reaching equilibrium. The slower merger shows that 2A/IT 
varies in large periods during the first moments of the evo¬ 
lution (in other words: it does not show a unique abrupt 
change in 2KIW, but rather two or more large periodic 
fluctuations before reaching equilibrium). Collapses, regard¬ 
less of being “cold” or “hot” show one large initial fluctu¬ 
ation amplitude in 2K/W. Interestingly, subsequent evolu¬ 
tion seem to be different: the “cold” collapse still experiences 
one more relatively significant fluctuation of 2A'/W before 
reaching equilibrium. The “hot” collapse, on the other hand, 
show a persistent, although of low amplitude, fluctuation of 
2K/W still for some time, when the “colder” collapse is 
comparatively well stabilized. 


3.5.2 The ‘Kandrup’ Effect 

In order to understand the dynamical behaviour presented 
by the simulated models, we apply a diagnostic advocated 
by Kandrup et al. (1993). The merging of stellar systems 
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occurs because of a transfer of the orbital energy to the par¬ 
ticles of the stellar systems in question. The mechanisms 
through which this occurs are the tidal interactions, which 
increase the internal energy of the systems at the expense 
of their orbital energy. The question here concerns the re¬ 
lation of this mechanism with the central non-homology of 
the simulated mergers. 

During the evolution of the system, the energy of the 
particles is not, in general, conserved, even in a “coarse 
grained” sense (viz. through the distribution N{E)-, for a 
discussion on th e importance of this distribution for stel¬ 
lar systems, see Binney 1982| ). Kandrup et al. studied the 
distribution of the energy of the particles in systems result¬ 
ing from collisions (without the formation of a final single 
object) and merging of two galaxies. These authors found 
that there is a “coarse-grained” sense in which the ordering 
of the mean energy of given collections of particles is un¬ 
altered, even though N{E) may vary substantially. In this 
section, we revisit the question raised by Kandrup et al. and 
try to connect this fact to the behaviour of the simulated 
systems in context of the FP. Notice that their conclusions 
were based on only two simulations of collisions, with only 
one merger, and two collapses. Here we use a much larger 
set of simulations and initial conditions, and a larger num¬ 
ber of particles as compared to the models used by Kandrup 
et al. We will not consider time evolutions of mean energies, 
as Kandrup et al. did, but only the initial and final values 
of the mean energy of the given collections of particles. We 
discuss further their diagnostic below. 

The method may be considered a “lagrangian” ap¬ 
proach to the analysis of how the energy of the particles 
change because of the relaxation process. The particles of the 
initial models have been sorted accordingly to their binding 
energies and the models were partitioned into 5 bins of equal 
number of particles (a finer partitioning with 10 bins pro¬ 
duced essentially the same results). For each of these bins, 
the mean energy was calculated and the bins ranked with the 
first one initially containing the most bound particles (most 
negative mean energy) whereas the fifth, the less bounded 
ones (less negative mean energy). The mean energy of these 
collections of particles were then recalculated at the end of 
the run’s and compared with their initial values. We have 
limited our analysis for the first generation of mergers. In 
the case of mergers with equal Eorb’s, we have included only 
the model with lower Lorb. The results of these comparison 
are displayed in Figs. and |^. 

We found that, except for some of the C cases, collapses 
preserve the ordering of the mean energies per bin entirely. 
These results confirm the findings of Kandrup et al. More¬ 
over the mean energies per bin that changed more in this 
case were the ones corresponding to the most bounded bins 
(1, 2, 3, etc). The central potential becomes deeper after the 
collapse, and the particles initially more bounded to the sys¬ 
tem tend to loose energy, becoming even more tied. This ef¬ 
fect is also a function of the collapse factor /3, as can be seen 
by the dashed line in the figure, which connects the most 
bounded bin (1), illustrating how this bin changes as a func¬ 
tion initial condition. In the case of A collapses, the mean 
energies representing the three less bounded bins converge to 
similar final values, whereas the two most strongly bounded 
bins reach even more negative mean values. This behaviour 
indicate that collapses tend to produce core-halo structures. 
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Figure 10. Energy means of 5 equal number collections of par¬ 
ticles, ordered accordingly to their initial mean binding energies. 
The initial and final mean binding energy values are connected 
by a line segment (left extreme: initial mean value; right extreme: 
final value), centered in the initial condition of each model. “Bin 
numbers” are indicated in the figure: bin number 1 refers to the 
most bounded particles; and successively throughout bin number 
5, which refers to the most weakly bound particles. The values of 
the final mean binding energy of the initially most bounded bin 
(1) are connected with a dashed line, illustrating how it changes 
as a function of the initial condition. The models shown are col¬ 
lapses. 

In the case of C collapses, the mean energies change con¬ 
siderably and chaotically. Recall that these models contain 
a Hubble flux which may favor the grow of perturbations 
embedded in these models, adding some complexity in the 
evolution of the mean energy of these distinct collections of 
particles. 

In the case of mergers (one and two-component mod¬ 
els), the preservation of the ordering of the mean energies 
per bin is not as good as in the case of collapses. For the 
D models (Hernquist one-component models), the initially 
most bounded particles will remain as the most bounded 
ones after the relaxation process. However, it can be seen 
that the bins nnmber 3 actually cross bins number 2 and 
reach the values corresponding to bins number 1. In the 
case of two-component models, the luminous matter tend 
to reach very negative values of the mean energies, almost 
converging to similar values for all bins. The general be¬ 
haviour of the luminous component resemble the behaviour 
of the most boundly tied particles of the collapse models. 
The main reason for this may be the fact that, after the 
initial interaction, the luminous component finds its equi¬ 
librium state within the deeper potential well of the dark 
halo. This might occur through a partial collapse of the lu¬ 
minous matter inside the dark halo. In the case of the dark 
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Figure 11. The same as previous figure, but for first genera¬ 
tion mergers. In the case of mergers with equal Eorb^^ we have 
included in the figure only the model with lower Lorb- 

matter component, there is also some violation of the order¬ 
ing of the mean energies per bin. In fact, the initially most 
bonndly tied bins (number 1) crosses upwards and gain en¬ 
ergy in some cases. The halo seems to be the only system 
that actually shows clearly this behaviour. 

The D models and the luminous component also present 
the same effect as seen in collapses: the particles initially 
more bounded to the system tend to loose energy as a func¬ 
tion of the initial orbital energy, as can be seen by the dashed 
line in the figure, which connects the most bounded bins (1). 
For halos, this behaviour is not as clear. 

Note that the mean biding energy of the most bounded 
bins remains at an almost constant value for collapses (see 
dashed lines in Fig. [lo| ), rising steeply for the “hottest” 
collapses. The C models present more fluctuations in this 
behaviour. For mergers, on the other hand, these changes 
proceed more smoothly and systematically with the initial 
orbital energy (dashed lines in Fig. pd] ). This means that 
the initially most bounded collection of particles remain at 
the average the most bounded particles after the relaxation, 
but at a more negative mean energy than the slower (less 
negative Eorb) mergers. 

It is clear that if the ordering of the mean energies of 
particles, partitioned at a coarse-grained level, is strictly 
conserved, as in collapses, then the most bounded particles 
(in average closer to the baricenter) continue be the most 
bounded particles after virialization. However, as pointed 
out previously, this ordering conservation does not occur for 
mergers. In other words, some complex behaviour seem to 
take place during the merger involving the more central or 
bounded particles, an effect which does not occur at all in 
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collapses (except for some small shufflinR between the less _ towards a core in less than a Hubble time (c.f. Yoshida et 


bounded bins for the “hottest” C models). In summary, we 


al 


2000). If some of these processes have operated in ha¬ 


entirely confirm the results of Kandrup et al. for collapses, 
but in the case of mergers, some violation of the ordering 
conservation of the energy bins is present. 


4 DISCUSSION 

We have simulated a hierarchical non-dissipative merger 
scheme similar to that of CdCC95, however using different 
models for representing the progenitors of the genera¬ 
tion mergers. In contrast with CdCC95, which considered 
King density profiles, we used models endowed with cuspy 
profiles, such as the Hernquist density profile. Also models 
with a dark halo second component were used in this study. 
A comparison with collapse simulations (previously anal¬ 
ysed in DCdCROl and here extended to include collapses 
with initial spins) is presented. 

We found that the one-component Hernquist mergers 
give results similar to those found by CdCC95 for the one- 
component King models, namely both were able to build-up 
small scattering FP-like correlations with slopes consistent 
to those found for the near infrared FP of nearby galaxies. 
The two-component models also reproduce a FP-like correla¬ 
tion, but with a signihcantly steeper slope which is in agree¬ 
ment with that found for galaxies at high redshift (Pahre 


1996). Pahre hnds that the slope of the near-infrared FP de¬ 


creases with increasing redshift (see his Figure 7.2). Another 
important piece of evidence of th e evolution of t he FP with 
redshift comes from the work of Kelson et al. (1997). The 


authors find that the structure of the galaxies in the anal¬ 
ysed sample has not changed significantly since a = 0.58, 
based on the fact that the observed scatter is rather low: 
±0.067 in log re. Besides, they Hnd a dependence between 
M/Lv and redshift, which reinforces the idea of a stellar 
population effect in the evolution of the FP. 

In nature, dissipational effects must have played a role 
in producing the FP relations and their scatter, but it would 
be unwise to completely disregard the role of stellar dynam¬ 
ics in shaping ellipticals as well. Our simulations are only 
of a dynamical character, with the M/L ratio fixed by con¬ 
struction. Systematic non-homology in the evolved models 
can produce FP-like ‘tilts’ compatible with those found in 
nature. In particular, we show the importance of the grav¬ 
itational potential of the halo for changing the ‘tilt’ of the 
FP: the magnitude of the change can be seen directly from 
a comparison between one-component and two-component 
Hernquist merger models (c.f. Fig. 1). In other words, this 
simple result clearly shows how the FP slope may be dynam¬ 
ically changed just by the addition of a halo. Therefore, our 
results suggest that the structure evolution of the halo could 
also have a collateral importance in changing and shaping 
the FP ‘tilt’, along with population evolution effects (e.g. 


Kelson et al. 1997). We speculate on the possibility that ha¬ 
los may suffer evolution from z ~ 0.5 to the present. The 
evolution could be in the form of tidal stripping, which would 
decrease the mass of the halo, or by the presence of super- 
massive central black holes, which could alter the matter dis¬ 


tribu tion of the halo, forming a core (c.f. Hennawi & Ostrikei 
2002 ). If dark matter presents some level of self-interaction 
(c.f. Spergel fc Steinhardt Mo|) , then it may drive evolution 


los during the past few Gyrs, in the sense of altering their 
gravitational potentials at the centers of galaxies, the char¬ 
acteristic scaling properties of the luminous/barionic com¬ 
ponent might have changed as well. Whether any of these 
possibilities are in factible is at present an open issue. 

A qualitative analysis of the behaviour of the mean en¬ 
ergy of collections of particles (as advocated by Kandrup 
et al. 1993) lead us to consider the possibility that ‘meso¬ 
scopic’ constraints could have some connection to the central 
non-homology. The conservation of the ordering of the mean 
energy of collections of particles implies that the process of 
‘mixing’ in the one-particle energy space is quite inefficient 
as compared with ‘mixing’ in conhguration and/or velocity 
space (see discussion in Kandrup et al.). This seem to be 
true for collapses, but not entirely for the central parts of 
mergers. 

The most intense tidal perturbations (shocks) seem to 
be found for the most rapidly merging systems (more neg¬ 
ative orbital energies). In this case, the particles probably 
withdraw the energy from the relative orbit of the merging 
pairs at one major fluctuation. Secondary fluctuations on the 
gravitational potential evolve afterwards rapidly and reach 
equilibrium in a short timescale as well. The stratification of 
the kinetic energies resemble that of the progenitor in this 
case. On the other hand, if the orbital energy is less nega¬ 
tive (slow mergers), there is some of time for the particles to 
withdraw energy from the orbit of the pair, and this process 
involves periodically large fluctuations on the potential that 
evolve slowly, taking a larger amount of time to stabilize. 
This process may be important in “heating up” the central 
parts of the models approaching Eorb —> 0. This should be 
important in defining the non-homology in mergers because 
the stratification of the kinetic energies are indeed different 
to that of the progenitor. 

On the other hand, in the case of collapses, the dynam¬ 
ics seem to operate in a different manner than in mergers. 
Collapses starts off from a spherically symmetric condition 
that mergers do not share. Collapses also produce not only 
fast by very high amplitude gravitational potential fluctu¬ 
ations that dump rapidly. This process should be very effi¬ 
cient in heating up the central parts of the models in con¬ 
figuration and/or velocity space, but not efficient enough to 
have the particles ‘forget’ their initial energies in a collective 
(‘mesoscopic’) manner. As already pointed out by Kandrup 
et al., this behaviour is at odds with Lynden-Bell’s theory 
of ‘violent relaxation’, where ‘mixing’ in energy space is not 
expected to be inefficient for any given collection of parti¬ 
cles. At the same time, we have found that collapses seem 
to ‘prefer’ forming homologous systems, whereas mergers 
do not. Some connection between ‘mesoscopic’ constraints 
and non-homology seem to be apparent, but this is an open 
question. 

We did not atempt at this time to rigorously try to 
connect the behaviour of the violation of the ordering of 
the energy bins (Figs. ^ and 0 ) with the behaviour of the 
‘central’ to ‘envelope’ kinetic energies (Figs. ^ and Al¬ 
though interesting, in order to fully understand this effect, 
we would need to probe the problem of relaxation in a much 
deeper and/or formal manner, which is not the objetive of 
our paper at the present time. Figs. ^ and ^ illustrate a 
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diagnostic on the behaviour of the change of energy of the 
CTstem due to relaxation in a ‘mesoscopic’ scale. Figs, ^and 
y show a different dignostic, where the change of energy of 
the central and external parts of the system are compared to 
their proqentor models, not the to their initial condition (as 
in Figs, and [l^) , and hence refer to a more ‘macroscopic’ 
feature of the relaxation process. On the other hand, the 
main point of Figs. and is not to show how the energy 
of collections of particles change due to relaxation (although 
it also certainly shows that) but in what degree the order¬ 
ing of their mean energies is violated. This type of analysis 
was first envisaged by Kandrup in 1993 and is still not well 
understood. In our oppinion, it is not at all clear how one 
could hud any immediate connection between both sets of 
figures. We know that the non-homology comes primarily 
as a systematic function of Eorb- Both sets of figures show 
systematic behaviour of 2 different types of diagnostics as 
a function of Eorb- Collapses show almost no dependecy of 
these same diagnostics with /3 (except for the “hottest” col¬ 
lapses). Therefore, it seems that relaxation through merging 
embodies some mechanism which is effective in differentiat¬ 
ing the final models, producing non-homology, whereas this 
mechanism is absent or highly precluded in collapses (again, 
except for the “hottest” ones, which are just a small pertur¬ 
bation form equilibrium of the progenitor model). In fact, 
what is lacking in order to make any progress in this di¬ 
rection is a systematic understanding of the nature of the 
gravitational relaxation mechani sm, where several co ncep- 
tual issues are still unsolved (c.f. Padmanabhan 1990). 

In any case, our results seems to strenghten the idea 
that dissipationless merging could produce signihcant non¬ 
homology in the final objects and therefore FP-like relations 
in the same sense and with comparable values of the FP 
‘tilt’ as those observed in ellipticals. We have shown that, 
from purely dynamical grounds, mergers can produce FP- 
like relations while simple collapses cannot (two-component 
collapses were not investigated here and will be a subject 
for future work). The evolution of gradients in the gravita¬ 
tional field of the merging galaxies seem to dictate the final 
non-homology of the end products. Further investigations 
are necessary in order to stablish, quantify and rigorously 
explain these complex effects, as preliminary discussed in 
this paper. 
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